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Abstract 

We present results for ground state structures of small Si n H (2<n <10) clusters 
using the Car-Parrinello molecular dynamics. In particular, we focus on how the 
addition of a hydrogen atom affects the ground state geometry, total energy and the 
first excited electronic level gap of an Si n cluster. We discuss the nature of bonding 
of hydrogen in these clusters. We find that hydrogen bonds with two silicon atoms 
only in Si2H, SisH and SisH clusters, while in other clusters (i.e. Su}i, SiqB., S17H, 
SisH, SigH and SiioH) hydrogen is bonded to only one silicon atom. Also in the case 
of a compact and closed silicon cluster hydrogen bonds to the cluster from outside. 
We find that the first excited electronic level gap of Si n and Si n H fluctuates as a 
function of size and this may provide a first principles basis for the short-range 
potential fluctuations in hydrogenated amorphous silicon. Our results show that 
the addition of a single hydrogen can cause large changes in the electronic structure 
of a silicon cluster, though the geometry is not much affected. Our calculation of 
the lowest energy fragmentation products of Si n H clusters shows that hydrogen is 
easily removed from Si n H clusters. 



PACS numbers: 71.15.Pd, 73.22.-f, 61.46,+w 



I INTRODUCTION 

During the last decade clusters have attracted a lot of attention because of their inter- 
esting and novel properties [ffj. Of particular interest are the binary clusters of hydrogen 
and silicon which are thought to be present in hydrogenated amorphous silicon(a-Si:H), 
porous silicon and silicon surfaces. In addition to the fundamental interest, their study 
may throw some light on complex phenomena occurring in these systems. Hydrogen plays 
an important role in these systems in phenomena like photoluminescence of porous sil- 
icon, potential fluctuations and the Staebler-Wronski effect in hydrogenated amorphous 
siHcon(a-Si:H)[|]-[16|. To understand these phenomena it is important to study how 



the addition of hydrogen affects the local electronic structure and geometry in these 
systems[@], [[§. Since these systems are very difficult to handle computationally, some 
understanding in this regard can be gained by simpler calculations on small hydrogenated 
silicon clusters. With this motivation, we have carried out a detailed study of ground 
state structures and electronic properties of small Si n H clusters( 2< n < 10) using the 
Car-Parrinello molecular dynamics(CPMD), focusing particularly on the effects caused by 
hydrogen. In our study we have investigated: (1) ground state geometries of Si n H clus- 
ters, (2) effect of hydrogen on the geometry of a silicon cluster, (3) stability of a silicon 
cluster due to the addition of hydrogen, (4) the first excited electronic level gaps of Si n H 
and Si n clusters, (5) bonding nature and position of hydrogen in silicon clusters and (6) 
the lowest energy fragmentation products of Si n H and Si n clusters. 

Several calculations have been done for many silicon hydrogen clusters by using various 
techniques. Using the Car-Parrinello method, Onida and AndreoniJlT]] studied the ground 



state geometry and electronic structure of hydrogen passivated crystalline fragments of 
silicon such as Si 5 H 12 , Si 6 H 16 , Si 8 H 8 , Sii H 16 and Sii 4 H 2 4. They found that Si-Si bond 
lengths were insensitive to size effects, but electronic properties were strongly affected. 
They also found that HOMO-LUMO gaps were not simply related to the size of the clusters 
and the localization of electronic states near the gaps was not necessarily silicon-like, even 
though the clusters are crystal fragments of silicon passivated by hydrogen. Quantum 



chemical calculations of Si n H m clusters were carried out by Meleshko et a/[|Tg] for n=6-16 
and m ranging from 2 to 20. They found that each H atom was bonded with only one 
silicon atom and localized outside the silicon skeleton and that the packing density in the 
skeleton decreased as the hydrogen content of the cluster increased. Miyazaki et a/[|H3[l 
performed density functional calculations for small hydrogenated silicon clusters of SigH^ 
(0< x <14) and showed that for the sequence Si6H :c _ 2 + H 2 — > Si6H x , the attachment 
of H occurred not at the site of silicon having dangling bonds but at the site where the 
LUMO of SieHa; has a large amplitude for x=2 and 6. According to this calculation, 
the bonding interaction of Is orbitals of hydrogen atoms with the LUMO of Si6H x _ 2 
should be the major cause of stabilization of the clusters. Their explanation makes it 
clear that hydrogen does not simply attach with silicon to saturate the dangling bonds, 
but it interacts at electronic level. 



Swihart et a/pU| used ah initio molecular orbital calculations to investigate struc- 



ture and energetics of selected hydrogenated silicon clusters containing six to ten silicon 
atoms. The clusters investigated were those that played the most important role in par- 
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ticle nucleationpTI] in silane during chemical vapor deposition. Shvartsburg et al [22 



modeled the dissociation of neutral and positively charged Si n clusters in n<26 range. 
They used dissociation energies to test the results of global optimization and fragmenta- 
tion products of the clusters. Recently, non-orthogonal tight-binding molecular dynam- 
ics(NTBMD) with simulated annealing optimization method was used to calculate ground 
state geometries of small Si„H clusters p3[| (2< n <10) and Si2H m (n=l, 2 and m=2-6) 



clusters [2^j. Using the non-orthogonal tight binding method, genetic algorithm opti- 
mizations were carried out for Si n H m clusters(n=l, 2 and m=2-6)p6], Experimental 
studies have been carried out for hydrogenated silicon clusters using a quadrupole ion 
trap[^8] where Si n H+(n=2-10 and x=0-20) were grown from silane gas. From the mass 
spectra of these clusters, it was shown how the stability of a silicon cluster is affected by 
hydrogenation. 

Our CPMD calculations show that hydrogen does not cause any drastic change in the 
geometry of the host silicon cluster although there is some distortion to the structure. 
To see clearly how the addition of a hydrogen atom affects the structure, stability and 
electronic properties of the host silicon cluster, we have also done a number of calculations 
for host silicon clusters. We have discussed two kinds of stabilities, one is geometrical 
stability and the other is electronic. To examine the geometrical stability of an Si n 



cluster |^9[-||32[, we have calculated the difference between the total energy of modified Si n 
geometry which is obtained by removal of hydrogen from Si n H cluster and ground state 
energy of Si n cluster. This energy difference gives information about the modification of 
the host silicon geometry due to the addition of hydrogen. To examine the electronic 
stability of Si n H cluster we have calculated the first excited electronic level gap for Si n H 
clusters. Comparison of the first excited electronic level gaps of Si„ clusters and Si n H 
clusters shows that hydrogen, in general, brings electronic stability to silicon clusters. We 
have also calculated the lowest energy fragmentation products of Si„ and Si„H clusters. 

The plan of the paper is as follows. In section II we give computational details of the 
present work. In section III the ground state geometries are presented and discussed in 
detail. In section IV we discuss the stability, cohesive energies, total energy differences 
between clusters, the first excited electronic level gaps and the lowest energy fragmentation 
products of Si n H and Si n clusters. Finally we summarize our results in section V. 



II COMPUTATIONAL DETAILS 

We have used the Car-Parrinello molecular dynamics(CPMD)P3|], pi with simulated 



annealing optimization technique to find the ground state structures of Si n and Si n H clus- 
ters. The CPMD method combines the density functional formalism with the molecular 
dynamics simulation. This scheme allows us to describe dynamics of ions under action of 
forces calculated by the Hellman-Feynman theorem. The pseudopotentials for silicon and 
hydrogen have been generated using the Bachelet, Hamann and Schl-uter technique p5 



The local density approximation(LDA) of the density functional theory has been used with 
the Ceperley- Alder [^] exchange-correlation energy functional parameterized by Perdew 
and Zunger [|37ll . 

The wave functions were expanded in a plane wave basis with 12 Rydberg energy 
cut-off and k=0 point was used for Brillouin zone sampling. During simulation volume 
of the system was kept constant and to avoid interaction between the clusters a big fee 
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supercell with side length of 35 a.u. was used. To perform simulated annealing, the 
system was taken to high temperatures (1200K in the steps of 300K), equilibiriated for 
a long time (about 16,000 steps) and then slowly cooled down (in the steps of 50K) to 
300K. Below this temperature the steepest descent optimization was found to be more 
efficient to obtain the ground state geometry. To check the upper limit of temperature, 
some of the clusters were heated up to 1500K and 2000K and it was found that the 
resulting ground state structures were the same. The desired temperature was achieved 
by rescaling atomic velocities and the atoms were moved according to the velocity Verlet 
algorithmic with a time step of 5 a.u. The fictitious mass of the electron was taken 
to be 200 a.u. All calculations were performed with more than one initial condition. 
The initial structures for MD calculations were chosen without any pre-assumption about 
the ground state geometry of the cluster. The starting atomic configurations were chosen 
arbitrarily with a constraint that atoms were neither very far away from each other nor too 
closePm , As mentioned above these clusters were then heated to high temperatures 



and then equilibiriated for a very long time. At this stage we find that the geometry of 
the hot cluster does not have any resemblance with the initial structure. At least two 
such CPMD calculations were performed for each cluster. For some clusters we have done 
CPMD calculation with three(SigH) and four( SigH) starting atomic configurations and 
found that the final structures are the same. Furthermore, we performed the CPMD with 
steepest descent optimizations on the NTBMD structures and found that the resulting 
geometries either converge to our structures or get stuck in some local minima. The 
ground state structures of Si n ( 2< n < 10) clusters obtained by us using the CPMD are 
similar to that obtained by previous calculations p9|, PD| . Also our result for bond length 



of SiH cluster(1.583A) is close to the earlier CPMD result^] and other calculations flO 
This shows that our calculational procedure is able to give correct structures. 

The first excited electronic level gap (e i+i-€{) of a cluster is calculated by transferring 
a small charge from its ground state configuration to its first excited state, and is given 
by[0 



e f+i-Ci = SE/ 5q 

where 5E is the difference between the total energy when 8q amount of charge is transferred 
to the first excited state and total energy of the ground state. 



Ill GROUND STATE STRUCTURES 



In this section we discuss in detail the ground state geometries of Si n H clusters obtained 



by the CPMD|33| with simulated annealing and steepest descent optimizations. By com- 
paring our ground state geometry of Si n H cluster with Si n cluster, we have investigated 
the effect of hydrogen on the geometry of the host silicon cluster. Furthermore, we have 
compared our ground state geometries of Si n H clusters with earlier calculated geometries 
of Si w F|4"2]| and Si n Na[^3| clusters. We have also made a detailed comparison of our work 
with earlier NTBMD |23| work. It agrees with our geometries of Si2H, SisH, Si4H and 



SiioH but the remaining geometries are different from our geometries, particularly the 



position and bonding of hydrogen. In the NTBMD |23| results, hydrogen was found to be 
bonded with more than one silicon in most of the clusters but in the present case we find 
this only for Si 2 H, Si 3 H and Si 5 H clusters. We find that our structures have lower energies 
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than those of the NTBMD structures. The nature of bonding has been investigated by 
performing charge density calculations. In the following we discuss our results for each 
cluster. 



(a)Si 2 H 

The ground state geometry of Si2H cluster is shown in Fig. 1(a). Two silicon atoms 
and hydrogen form a triangular structure. Note that the two silicon atoms are bonded 
to each other not only via Si-Si bond but also via Si-H-Si bridge type bond. The lowest 
energy structure of Si 2 F||2|] and Si 2 Na||3] are similar to the structure of Si 2 H. In Si 2 H 



cluster, Si-Si bond length is 2.131A and hydrogen is equidistant from both silicon atoms 
with the bond length of 1.724A which is larger than its value of 1.583 A in Si-H dimer. 
The Si-Si bond length in Si 2 H cluster is smaller than the Si 2 dimer bond length of 2.184A. 
This implies that hydrogen pulls both silicon atoms closer and increases bonding between 
them. Thus, the additional bonding between silicon atoms is due to the Si-H-Si bridge 
bond which is attractive in nature. Such bridge type bonds are thought to be present 
in a: Si-H and play an important role in explaining the Staebler-Wronski effect H , [BJ , |T2 
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It is interesting to note that hydrogen is bonded with both silicon atoms although 
its valance is one. This is seen clearly from the valance charge density plotted in Fig. 
2, which shows existence of bonds between hydrogen and the two silicon atoms. Such 



overcoordination of hydrogen has also been observed recently in SiC system f4"5| . We also 
see from this figure that the electrons are more localized near the hydrogen atom. This 
is expected since hydrogen is more electronegative than silicon Thus there is a small 
charge transfer from silicon atoms to hydrogen |TJ] . As a result, Si-H bond is neither 
purely covalent nor ionic but is polar covalent p8| . 

(b)Si 3 H 

The ground state geometry of SiaH is shown in Fig. 1(b). This is a planar structure 
with two fold symmetry and has some resemblance with Si4 cluster 0]. Hydrogen in this 
cluster is bonded with two silicon atoms (1 and 2) and is equidistant from both atoms 
with bond lengths of 1.715A. Silicon atom No. 3 is also equidistant from silicon atoms 1 
and 2. The bond length between silicon atoms 1 and 2 is 2.377A which is more than 
the Si-Si bond length in Si 2 H. This indicates that bonding between 1 and 2 silicon atoms 
is weaker than the Si-Si bond in Si 2 H because of the presence of another silicon atom. 
Comparing this with Si 2 H structure, we note that the additional silicon takes diagonal 
position opposite to hydrogen. We see that although hydrogen does not modify the Si 3 
geometry much, it does modify the bond lengths. Particularly, the bond length between 
silicon atoms 1 and 2 in SisH is smaller than bond length of 2.613A in Si3 cluster. Thus 
as in Si 2 H, hydrogen pulls silicon atoms 1 and 2 closer, which can be attributed to the 



Si-H-Si bridge bond. We note that the lowest energy structure of Si3Na|43[ is similar to 



the present structure, but for Si 3 F cluster the ground state geometry is different |4"2"| . 
(c)Si 4 H 

The ground state geometry of the cluster is shown in Fig. 1(c). Four silicon atoms 
form a flat rhombus and the hydrogen atom is above the plane and bonded with one 
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of the silicon atoms. The same structure was shown as the lowest energy geometry of 
Si 4 F cluster The lowest energy structure of Si 4 Na|43| cluster is similar to the present 



structure but differs in the coordination of Na atom. Comparison of this structure with 
Si 4 structure pOf shows that the addition of hydrogen does not bring much change to Si 4 
structure. Comparison with Sis cluster p0|] shows that this structure does not have any 
resemblance with Si 5 cluster. Based on the idea of local softness and hardness, Galvdn et 
q/f49| predicted the sites preferred by hard and soft species in Si 4 cluster. We find that 
hydrogen goes to the position according to their prediction. 

(d)Si 5 H 

Two lowest energy structures of SisH are very close in energy and differ only by 0.06 
eV. The geometry which is higher in energy is shown in Fig. 1(d). Three of the silicon 
atoms, numbered 1, 4 and 5 form a triangular plane and silicon atom No. 2 is above 
and No. 3 is below the plane. Hydrogen takes the apex position in the structure and is 
bonded with only one silicon atom. The geometry of the silicon atoms is same as in Sis []30f 



cluster. The lowest energy structure of SisH cluster is shown in Fig. 1(e). Geometry of 



this cluster is similar to Si6 cluster |30| except that one of the silicon atoms is replaced by 
hydrogen. Note that hydrogen in this cluster is attached to two silicon atoms, which is 
also the case in Si2H and SisH clusters. Geometrically hydrogen plays the role of silicon 
in these three clusters, i.e. the geometry of Si 2 H is similar to Si 3 , Si 3 H is similar to Si 4 
and SisH is similar to Si$. We note that a two coordinated silicon atom exists in Si3, 
Si 4 and Si 6 clusters and hydrogen replaces this silicon atom to form Si 2 H, Si 3 H and Si 5 H 
respectively. We speculate that this may be a general feature of Si n H clusters, i.e., if 
a two coordinated silicon atoms exists in a Si„ + i cluster, hydrogen will replace the two 
coordinated silicon atom to form Si n H cluster, which will have the same geometry as 
Si n+ i cluster. This seems to imply that hydrogen will form Si-H-Si bridge type bond 
between two nearby silicon atoms which are doubly coordinated. In the context of a-Si:H 
this would imply that hydrogen will form Si-H-Si bond between two nearby silicon atoms 
having two dangling bonds. Comparison with SisF geometry shows that one of the low 



energy structures j42] is similar to SisH shown in Fig. 1(d). The low energy geometries for 



SisNa[43"] cluster are similar to our low energy geometries for SisH and the lowest energy 
structure is also the same. 

(e)Si 6 H 

The ground state geometry of Si 6 H is shown in Fig. 1(f). In this structure, four 
silicon atoms numbered 1, 2, 3 and 5 form a distorted plane and the remaining two atoms 
numbered 4 and 6 are above the plane. Hydrogen atom is bonded to silicon No. 4 from 



outside the cluster. Comparing this with the ground state geometry of Sie [^0[ cluster, we 
note that although the plane formed by four silicon atoms is same as in Sig cluster, the 
other two silicon positions are different i.e. in Sig[|3U] cluster one silicon atom is above 



and another is below the Si 4 plane, but in SigH two silicon atoms are above the plane. 
Comparing with the results of SiqRx clusters jn|, our geometry of Si^H. falls in the class 
of tetrahedral bonding network. We find that this is the only cluster among the clusters 
considered here where the geometry of the silicon atoms differs from the ground state 
geometry of the host silicon cluster pOj], This shows that hydrogen can cause a transition 



from one geometry to another geometry. One of the low energy geometries of SigF| 
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cluster is same as the present structure. But in the case of SigNa[|^] one of the low energy 
structures has similar geometry but it differs from our structure by coordination of the 
Na atom. 



(f)Si 7 H 

The ground state geometry of the structure is shown in Fig. 1(g). Silicon atoms in this 
structure form a closed and compact unit and the cluster has pentagonal symmetry. Five 
silicon atoms numbered 1, 7, 6, 4 and 3 make a pentagonal plane and one silicon atom 
is above and another is below the plane. Silicon atoms which are not in the pentagonal 
plane are bonded to all the atoms in the pentagonal plane. Hydrogen takes the apex 
position in the structure and is bonded with one silicon which is out of the pentagonal 
plane. It is interesting to note that instead of bonding with four coordinated silicon, 
hydrogen is bonded with five coordinated silicon. In Si 7 H, Si 8 H, Si 9 H and Sii H we found 
the same trend of hydrogen preferring to bond with an over-coordinated silicon atom. 
This is surprising since one would have expected it to bond with less coordinated silicon. 
This may be attributed to slightly higher electronegativity of H compared to Si and as 
a result, H prefers to bond with silicon atom having more number of electrons. This is 
consistent with the earlier calculation on structural evolution of SigHj. clusters [[0|], where 
it was found that hydrogen is not necessarily bonded with a silicon site having dangling 
bonds, but with a site where LUMO amplitude is larger. Comparison of Si 7 H with Si 7 



| 29|j , f30|j cluster shows that hydrogen hardly changes the geometry of Si 7 cluster implying 
that Si 7 cluster is a very stable cluster. The lowest energy geometry of Si 7 F[|42| and one 
of the low energy structure of Si 7 Naj|3] is similar to our ground state geometry. 



(g)Si 8 H 

Figs. 1(h) shows the ground state geometry of SigH. This structure is also made of a 
compact, closed unit of silicon atoms with hydrogen sticking to the structure from outside. 
We see that SigH shows some similarity with Si 7 H cluster i. e. the same pentagonal plane 
formed by silicon atoms numbered 4, 3, 1, 7 and 2 exists in Si 7 H and one silicon is above 
and another is below the plane(6 and 5) as in Si 7 H. Silicon atom 8 is attached to the 
triangular plane of the Si 7 H structure in such a way that it is away from hydrogen atom. 
As in Si 7 H hydrogen atom is attached to the silicon atom which is bonded with five silicon 
atoms. Fig. 1(h) shows that SisH has two distorted Si4 planes; silicon atoms numbered 8, 
4, 3 and 5 form one distorted Si4 plane and 2, 6, 1 and 7 form another. These planes are 
not parallel but rotated with respect to each other in such a way as to have more than one 
bond for each silicon atom with the atoms in the other plane. Comparing our ground state 
geometry of SigH with Sig|29[] geometry, we see that silicon atoms have similar geometry 



except that the Si4 planes are distorted in SigH. 
(h)Si 9 H 

The ground state geometry of the cluster is shown in Fig. l(i). This structure is also 
compact and closed by silicon atoms. The structure consists of two Si4 planes formed by 
atoms 1, 3, 4, 9 and 5, 8, 2, 7 and silicon No. 6 forms a cap. Hydrogen is connected to 
silicon atom No. 9 which is coordinated with five silicon atoms. Comparison with Sig[[2S 
cluster shows that the atoms forming Si4 planes in SigH do not lie in a plane in Sig. 
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(i)SiioH 



The ground state geometry of this cluster is shown in Fig. l(j). The Siio structural 
unit in SiioH is similar to SigH cluster and the additional silicon(atom 1) makes a side 
cap to SigH cluster ( Fig. l(j)). Also there are two Si4 planes rotated with respect to each 
other as in SisH and SigH. Silicon atoms 5, 2, 7 and 6 form one plane and 8, 3, 10 and 4 
form other plane. In this cluster hydrogen is connected to the five fold coordinated silicon. 
Comparison with Siio [^9| cluster shows that hydrogen hardly changes the geometry of 
Siio cluster implying that Siio cluster is a very stable cluster. A general feature of clusters 
SieH, SirH, SisH, SigH and Si^H is that silicon atoms in the cluster form a closed compact 
unit with hydrogen outside this structural unit. 



IV STABILITY OF Si n AND Si n H CLUSTERS 



We find that the total energy of Si n as well as of Si n H clusters increases approximately 
linearly with the cluster size n. The addition of hydrogen to an Si n cluster reduces the 
energy of the cluster by approximately 15 eV. The cohesive energy per particle versus 
number of silicon atoms is plotted in Fig. 3. As seen clearly from the figure cohesive 
energy per particle increases rapidly up to SieH cluster and then it increases slowly as 
a function of size. As noted earlier, from this size ( SieH) onwards silicon atoms in the 
cluster form closed compact unit and some of the silicon atoms have coordination number 
more than four. 

We take the first excited electronic level gap of a cluster as the difference between 
the first excited electronic level and the highest occupied level. For closed-shell or sub- 
shell systems this gap will be the same as HOMO-LUMO gap which is related to the 
chemical hardness and electronic stability of a system 1 50] -[^2]. Though the first excited 



electronic level gap is not equivalent to HOMO-LUMO gap for Si n H clusters, it can be 
related to electronic stability. A bigger value of the first excited electronic level gap for 
a system means that it is difficult to excite electrons from its ground state and thus the 
electronic system can sustain its ground state for larger perturbations. Thus the first 
excited electronic level gap can be taken as a measure of the electronic stability of a 
system. We have shown the first excited level g function of cluster size in Fig. 4 



for Si n H and Si n clusters. Also shown in the figure are results of Lu et al|30| for Si n clusters 
which are in good agreement with our results. We see that general trend of variation of 
first excited electronic level gap is quite similar for Si n H and Si n clusters. The figure also 
shows that the addition of hydrogen can cause large changes in the electronic structure 
of Si n cluster. 

From Fig. 4 we see that the gap fluctuates with size, which indicates that the gap 
strongly depends on the size and geometry of a cluster. It might be interesting to draw 
parallels with short-range potential fluctuations in a-Si:H system which occur at the length 
length scale of 3A[]I6 . It can be argued that an amorphous system can be considered 



as a loosely connected network of small clusters and thus our calculation provoids a 
first-principle basis for the potential fluctuations Oj-|TS]. Furthermore, we see that 



the first excited electronic level gap for Si n H is, on an average, larger than that of Si n 
cluster. This is consistent with the observation that the band gap of a:Si increases on 
hydrogenationp3| . Further, Fig. 4 shows that Si 2 H, Si 3 H, Si 5 H, SiyH, Si 9 H and Sii H 



clusters are electronically more stable compared to Si 4 H, Si 6 H and Si 8 H clusters. Also we 
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see that Si2, Sis, Sie, S17 and Siio clusters are electronically more stable than other silicon 
clusters(Si 3 , Si 4 , Si 8 and Si 9 clusters), since they have larger gaps. 

To examine geometrical stability we have calculated the difference between the total 
energy of the modified Si n cluster, which has the same positions of silicon atoms as in the 
Si„H cluster, and the ground state geometry of Si n cluster. This energy difference is a 
measure of how much a silicon cluster distorts from its ground state geometry due to the 
addition of hydrogen atom. Lower value of this difference for a Si n cluster means that the 
cluster is geometrically stable. This total energy difference as a function of cluster size n 
is shown in Fig. 5. The figure shows that Si 2 , Si 4 , Si 7 and Siio clusters are geometrically 
more stable than Si 3 , Si 5 , Si 6 , Si 8 and Si 9 clusters. On the other hand Si 3 , Si 5 , Si 6 , Si 8 and 
Si 9 clusters are stabilized by hydrogen and have a greater tendency to adsorb hydrogen. 
This is consistent with the conclusions drawn by comparing Si n H and Si„ ground state 
geometries in section III, as Si 2 , Si4, S17 and Sii clusters were least distorted by the 
addition of hydrogen. Thus the above discussion shows that, Si 2 , Siy and Sii clusters are 
the most stable clusters from both viewpoints of electronic as well as geometrical stability. 

As pointed out earlier in section III, hydrogen is attached to silicon clusters from 
outside in several cases (Si4H to SiioH). To examine this further, we have performed 
the steepest descent optimization on SieH and Si7H clusters with hydrogen surrounded 
by silicon atoms. We find that hydrogen atom always comes out of the silicon cluster 
independent of the cluster size. This is mainly due to the higher electrostatic energy of 
the cluster when hydrogen is inside the cluster. Thus our result implies that hydrogen 
will tend to come out of crystalline silicon and would like to stay on the surface. This 
is consistent with the experimental observation in which hydrogen is used to produce 
homogeneous silicon surface by terminating surface silicon dangling bonds to reduce the 
surface reconstruction fT!2|j . 

To investigate fragmentation products of Si n and Si n H clusters, we have calculated 
the difference between the total energy of a cluster which undergoes fragmentation and 
that of its possible product clusters. The most probable pathway for fragmentation of 
a particular cluster is the one which has the smallest total energy difference [^TJ. Since 
the clusters are small in size, we are assuming that the fragmentation results in only two 
product clusters. Our calculations are only for neutral fragmentation of Si n and Si n H 
clusters. Our lowest energy fragmentation products of Si n clusters agree very well with all 



primary fragmentation products calculated by Shvartsburg et a/[^2|. In Table. 1 we have 



given the lowest energy fragmentation products of Si n H clusters with the corresponding 
dissociation energies. We see from the table that the lowest energy fragmentation products 
have hydrogen atom as one of the products for all Si„H clusters except for Si 8 H cluster. 
This shows that it is easy to remove hydrogen from Si„H clusters. 



V CONCLUSIONS 

We have presented detailed results for the ground state structures and electronic properties 
of Si n H clusters using the Car-Parrinello molecular dynamics simulations. We find that 
hydrogen can form a bridge like Si-H-Si bond connecting two silicon atoms. Such bridge 
like bonds are thought to be present in a-Si:H[EJ|, [ pHfl . However, among the clusters 
considered here hydrogen forms a bridge like bond only in Si 2 H, Si 3 H and SisH clusters; 
in others it is bonded with only one silicon atom and attached to the cluster from outside. 
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Charge density calculations show that the Si-H bond in all clusters is polar covalent. In 
clusters from Si 7 H to SiioH, silicon atoms form a compact unit and hydrogen attaches 
to a silicon atom which is over-coordinated. Though hydrogen has small effect on the 
geometry of the host silicon cluster, it changes bond lengths and tries to distort the 
silicon cluster. This is similar to the behavior of hydrogen in a-Si:H where it has been 
found that hydrogen creates local distortions as it moves 0, [[|. We find that hydrogen 
has a tendency to come out of compact silicon clusters and prefers to stay out of the 



cluster. This is consistent with the behavior of hydrogen on silicon surfaces] 12 



The first excitation electronic level gap of the Si n H clusters fluctuates as a function of 
size and this may provide a first principles basis for the short range potential fluctuations 



in a-Si:H|T3[-JT6|. Our calculations show that the addition of hydrogen can cause large 
changes in the electronic structure of host Si n cluster. Furthermore, it shows that Si2H, 
Si 3 H, Si 5 H, Si 9 H and Sii H clusters are electronically more stable than Si 4 H, Si 6 H 

and Si 8 H clusters. We find that Si 2 , Si 4 , Si 7 and Siio clusters are geometrically more 
stable than Si 3 , Si 5 , Si 6 , Si 8 and Si 9 clusters, while Si 2 , Si 5 , Si 6 , Si 7 and Sii clusters are 
electronically more stable than Si 3 , Si 4 , Si 8 and Si 9 clusters. We have calculated the lowest 
energy fragmentation products of Si n and Si n H clusters. Our results for fragmentation 
products of Si„ clusters agree very well with the earlier predictions. The lowest energy 
fragmentation products of Si„H clusters show that it is easy to remove hydrogen from 
silicon clusters. Comparison of Si2H to S17H. cluster with the corresponding F and Na 
substituted clusters shows that almost all have similar low energy geometries implying 
that the geometrical effect of H, F and Na on Si n clusters are similar. 
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Table. 1. Fragmentation pathways of a neutral Si„H cluster into products Si m and 
SL_ m H cluster. 
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FIGURE CAPTIONS 



Fig. 1. Ground state geometry of (a) Si2H, (b) SisH, (c) Si4H. (d) Higher energy 
geometry of SisH. (e) Ground state geometry of SisH, (f) SieH, (g) S17H., (h) SisH, (i) 
SigH and (j) SiioH cluster. Silicon atoms are numbered and hydrogen atom is shown by 
a small dark circle. A thick line between two atoms indicates a bond between the atoms. 

Fig. 2. Valence charge density of Si2H cluster in arbitrary units plotted in the plane 
of the cluster. Constant charge density contours are also shown. Approximate positions 
of Si and H atoms are indicated by arrows. 

Fig. 3. Cohesive energy per atom of Si n H cluster versus number of silicon atoms. 

Fig. 4. First excited state electronic level gap of Si n H and Si n cluster versus number 
of silicon atoms in the cluster. The circles and squares correspond to Si n H and Si„ clusters 
respectively. The triangles represent the results of Lu et alffi^ f° r Si n cluster. 

Fig. 5. Difference dE, between total energy of the modified geometry of Si n and the 
ground state energy of Si n cluster versus number of silicon atoms in the clusters. 
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